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Abstract 

cS ■ 

We investigate in this article the long-time behaviour of the solutions to the energy- 
dependant, spatially-homogeneous, inelastic Boltzmann equation for hard spheres. This 
model describes a diluted gas composed of hard spheres under statistical description, that 
dissipates energy during collisions. We assume that the gas is "anomalous" , in the sense 
that energy dissipation increases when temperature decreases. This allows the gas to 
cool down in finite time. We study existence and uniqueness of blow up profiles for this 
model, together with the trend to equilibrium and the cooling law associated, generalizing 
i-^h , the classical Haff's Law for granular gases. To this end, we investigate the asymptotic 

behaviour of the inelastic Boltzmann equation with and without drift term by introducing 
new strongly "nonlinear" self-similar variables. 

'■ 1 Introduction 

> 

£T) ■ We are interested in this paper in a particular model of granular gases. A granular gas is 

a set of particles which interact by energy dissipating contact interactions. This is a quite 
different model than perfect molecular gas where energy is conserved, because the particles are 
"microscopic" regarding the scale of the system, but "macroscopic" in the sense that there are 
not molecules at all, but rather grains of a given size. For example, a suspension of pollen in 
a fluid or a planetary ring |21|, [3] can be seen as a granular gas when intermolecular forces are 
neglected, each one in a different scale. 

Our goal is to investigate the global behaviour of the kinetic energy of a space-homogeneous 
gas of inelastic hard spheres interacting via binary collisions (that is a granular gas), and 
existence and uniqueness of blow up profiles associated. The study of this kind of system 
started with the physics paper |2Uj of P.K. Haff, and has since generated a large increase of 
interest, both in Physics (a complete introduction of the subject can be found in the textbook 
|12] by N. Brilliantov and T. Poschel) and Mathematics (cf. the review article |3Uj by C. 
Villani). It involves some complicated phenomena, such as kinetic collapse (cf. K. Shida |28J) 
where the gas freezes completely, clustering at hydrodynamic level for inhomogeneous gases, or 
even spontaneous loss of homogeneity (as proved by I. Goldhirsch and G. Zanetti in [E]). This 
article will especially deal with the case of a gas whose particle's collision rate increases with 
dissipation of energy. Such a gas is sometimes called anomalous (see article |26| of S. Mischler 
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Figure 1.1: Geometry of collision (dashed lines are elastic and solid lines inelastic) 

and C Mouhot), because of the unusual behaviour of this collision rate. Another model of 
anomalous gas has been introduced by G. Toscani in [29] by assuming that the collisions are 
close to be elastic for large relative velocities of the grains. It has been studied precisely by H. 
Li and G. Toscani in [23J, showing the finite time extinction of the support. 

A granular gas can be described in a purely Newtonian way, but the number of macroscopic 
particles involved (ranging from 10 6 to 10 10 ) leads to adopt a statistical physics' point of view. 
Therefore, we shall study the so-called space-homogeneous inelastic Boltzmann equation (also 
known as granular gases equation), given by 

(1.1) % = Qe(f,f), 

where / = fit, v) represents the particle's distribution function, depending on time t > and 
velocity v G R . It is always a nonnegative function. The collision operator Q e (f, g), which will 
be made precise in the following (particularly concerning the meaning of subscript e), models 
a binary inelastic collision process of hard spheres type, localized in time. This equation will 
be supplemented with the initial value 

(1-2) /(0,-)=/m- 

The inelasticity is characterized by a collision mechanics where mass and momentum are 
conserved and kinetic energy is dissipated. At microscopic level, it can be described as fol- 
lows: given two particles of pre-collisional velocities v and v*, their respective post-collisional 
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velocities, denoted by v' and are given by 

V =V — {U-CJ)LU, 

/ . 1+e . , 
= v* -\ — (u ■ U))U, 

where u := v — is the relative velocity of a pair of particles, u is the impact direction and 
e 6 [0, 1] the dissipation parameter, known as restitution coefficient (see Figure [1] for a sketch 
of the collision process). 

Physically, it means that energy is dissipated in the impact direction only. The parameter e 
can depend on relative velocity and kinetic energy of the particles. For the need of the analysis, 
we will simply assume that e is a constant, which is a rather good physical approximation (a 
complete discussion on this topic can be found in |12j). Besides, e will also be taken different 
from 1, since e = 1 concerns the classical elastic case, in which no dissipation occurs. The 
inelastic case e < 1 can also be characterized by the fact that the impact direction oj does not 
bisect the angle between pre and post-collisional velocities. 

The parametrization of post-collisional velocities can also be found by using some properties 
of the model. Indeed, it is equivalent to conservation of impulsion and dissipation of energy: 

v' + = v + v*, 
1 - e 2 

\v'\ 2 + Kl 2 - \v\ 2 - K| 2 = — |u • oj\ 2 < 0. 

One of our goal will be the study of the macroscopic dissipation of energy, sometimes 
referred to as the cooling process of the gas. This is of high interest when one wants to study 
the self-similar behaviour of dissipative kinetic equations. 

This phenomenon has been investigated mathematically by S. Mischler and C. Mouhot 
in [24j for the case of constant restitution coefficient e and by R. Alonso and B. Lods in [2] 
for visco-elastics hard spheres, more realistic at the physical level. In the latter one, energy 
dissipation is a power law of the relative velocity of colliding particles. We will use ideas of these 
two articles to write some of the results we present in the following. The paper [23] also deals 
with the cooling process of a one dimensional nonlinear friction equation, which arises when 
dealing with quasi-elastic limit of equation (|l.ip with variable restitution coefficient. Finally, 
J. A. Carrillo, M. Di Francesco and G. Toscani in [13] and J. A. Carrillo and J.L. Vazquez in [2] 
used the behaviour of the kinetic energy to investigate intermediate asymptotics of nonlinear 
diffusion equation. 

Let us now describe precisely the collision operator we consider throughout the rest of this 
article. Let / and g be two nonnegative particles distribution functions only depending on 
v 6 M. d . The collision operator Q e (f,g), where e € [0, 1) is the constant restitution coefficient, 
can be expressed in the following weak form: given a regular test function tp, 

(1.3) (Q e (f,g),i>) ■= \ \ f*9 (V^ + V'* -ip-ip*) B(\u\, u-uj,£{f))du}dv dv*, 

z JR d xR d xS d ~ 1 

where we have used the usual shorthand notation tp' :=ip(v'), ip* :=tp(v'^), u := u/\u\ and 

£(/):= / \v\ 2 f(v)dv 
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denotes the kinetic energy of /. Moreover, B is a positive function known as the collision 
kernel. We will assume that it factorizes as: 

(1.4) B(\u\,u-u),£) = \u\b(u-u>)£- a , 

for a nonnegative constant a and a nonnegative function b (the collisional cross section) of 
mass 1 in the unit sphere, bounded by below and above by two positive constants fi\ and fa- 

(1.5) Vx G [-1, 1], < /3i < b{x) < fa < oo. 

The gas is anomalous thanks to this collision kernel. Indeed, according to (|1.4|) . the collision 
frequency will increase when the kinetic energy decreases. As we will see in the following, this 
kinetic effect will leads to the cooling in finite time, provided that the coefficient a is big 
enough. Such a phenomenon is related to the one described in [29} 123] : the gas is anomalous 
because the grains are close to be elastic for large relative velocities. This is in opposition with 
the more physically intuitive models where collisions are elastic for small relative velocities 
(for example viscoelastic hard spheres, as described in [12\ [2] or even the simplified threshold 
model of T. Poschel, N.V. Brilliantov and T. Schwager of [27]). It also yields a finite time 
cooling. 

Concerning this topic, we can also mention some recent works of I. Fouxon et al. |17| and 
Kolvin et al. |22j . Under a weak inelasticity assumption with constant restitution coefficient, 
these authors derived granular hydrodynamics (of Euler and Navier-Stokes type) and proved 
the finite time cooling in these models, for a one parameter family of initial conditions. To do 
so, they used a Lagrangian formulation to exhibit non self-similar solution to these equations, 
and we shall not consider this approach in our paper. 

We can also give a strong form of the collision operator. As it has been pointed out in 
|30| (the full derivation of this expression can be found e.g. in articles [5j[6] of A. Bobylev, J. 
Carrillo and I. Gamba), one can write 



Q e (f, g)(v) = 8(f)- a / \u\ 7 'g* - fg* b(u ■ a) da dv* 

JWd x §d-l \ \U\ J 

= Qt(f,g)(v)-f(v)L(g)(v). 

In this expression, 'v and 'u* are the pre-collisional velocities of two particles of given velocities 
v and u*, defined for a G S^ -1 as 

. v + v* 1 — e . . 1 + e . . 

v = — — (v - n*J H — \v - cr 

2 4e 4e 

. v + v* 1-e, , 1 + e. . 

v * = — 1 : — {v ~ v *) : — \ v ~ v *\°~- 

2 4e 4e 

The parameter a is no longer the impact direction, but the center of what we can call the 
collision sphere [30]. The operator Qt(fi9)( v ) is usually known as gain term because it can 
be seen as the number of particles of velocity v created by collisions of particles of pre-collisional 
velocities 'v and whereas f(v) L(g)(v) is the loss term, modeling the loss of particles of pre- 
collisional velocities v. Finally, the coefficient J is the Jacobian of the map (v,v*) — > ('v, 'v*), 
not identically equal to 1, but since e is a nonzero constant given by 
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Taking successively ip( v ) = 1) v an d \v\ 2 in (|1 -3j) gives macroscopic conservation of mass 
and momentum, and dissipation of kinetic energy. This implies that equilibria of this collision 
operator are Dirac distributions p5 v=vo of prescribed mass p and momentum pvQ. In order to 
prevent apparition of these trivial solutions, one can add a diffusive term to give an input of 
energy, as was studied for example by Cercignani, R. Illner and C. Stoica in [15] and by A. 
Bobylev and Cercignani in [7] for Maxwell molecules. It is also possible to look for self-similar 
solutions to this equation thanks to a rescaling which leads to study the inelastic Boltzmann 
equation with drift term, as was done by M. Ernst and R. Brito in |16| . 

Let us now describe more precisely the asymptotic behaviour of the energy of /, assuming 
that this distribution is time dependant, of mass 1 and zero momentum. 

Definition 1. If / is solution to (jl.ip . we call cooling process the asymptotic behaviour of 
£{f){t). We say that there is a blow up when £{f){t) — > for t — > T c , where T c is the time of 
explosion of /. 

If one multiplies equation (jl.ip by \v | 2 and integrates in velocity, using the weak formulation 
(jl.3|) of the collision operator and (\1A\) , one gets dissipation of kinetic energy: 



(1.6) 



d_ 

dt 



£ (/)(*) = -£{mr a 



f(t,v)f(t,v*)\u\D (|u|) dvdv*, 



where D is a nonnegative quantity usually called dissipation rate, given by 



D(\u\ 



1 



\u ■ ui\ 2 b(u ■ to) du. 



Using the polar coordinates cos(#) =u-lo, we have 



(1.7) 



D(\u\) 



1 



u ■ 


u, 


we 


S c 


1-2 


/ 






J 



cos\9)sm d - 6 {e)b{cos{6))d6. 



Provided that / is of mass 1 and zero momentum, one can then apply Jensen's inequality 
to ^{v) = |u| 3 and to the probability measure f(v^)dv^, together with conservative properties 
of equation (jl.ip to get 



/ f(v*)V(\v- v*\)dv* > * 
Using equalities 



v / f(v*)dv* 



/ w*/(f*)di 

JK. d 



*(!«!)■ 



(1.8) 

where /3 e is given by 



and (|1.7|) . it comes that 
d 



dt 



£{f)(t)<-p e 8{f)(t)- a / f(t,v)\v\"dv, 



§ d-2 


I 







cos 2 {e)sm d - 3 (e)b(cos{e))de < +oo. 



Moreover, applying Holder inequality with s = 3/2 and its conjugate s' = 3 to the map 
v i — y f(v)\v\ 2 , one gets thanks to mass conservation 



(1.9) 



f(v)\v\ 2 dv 



< 



f(v)\v \ 3 dv. 
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Using Fubini Theorem for nonnegative measurable function and gathering inequalities (|1 .8|) 
and (|1.9p . we finally find a closed differential inequality for £ (/), namely 

(1.10) Jt £U){t) ~ -^S(f)(tr a+ ^, for t < T c . 

According to this inequality, the cooling process will depend on the nonnegative coefficient 
a introduced in the definition of the collision kernel (|1.4p . We now have to describe rigorously 
the spaces of solutions to (jl.ip . in order to introduce the corresponding Cauchy problem and 
the main results of this article. 



1.1 Functional Framework and Main Results 

Let us present some functional spaces needed in the paper. We denote by L q for q G [l,+oo) 
the following weighted Lebesgue space 

L\ = If : R d -> R measurable; \\f\\ L i := j \f{v)\ (v) q dv < ooj , 

where (v) := \J\ + |v| 2 . Thanks to this definition, we can introduce the sets of distribution of 
given mass 1, zero momentum and energy £ > as 

G:=\f€L\: I f(v)dv = l,[ f(v)vdv = o) 

and 

Qe :={f£G:J^f(v)\v\ 2 dv = £y 

Then, we denote by W g ' the weighted Sobolev space 

W 1 / :={/€Lj;V/€Lj}. 

We also introduce the space BV q of weighted bounded variations functions, given as the set 
of weak limits in V of sequences of smooth functions which are bounded in Wq' 1 . Finally, 
A / i 1 (M (i ) is the space of probabilities measures on M. d . 

The Cauchy problem (jl.ltil.2p has been extensively investigated in the past few years, 
starting with paper [5] and then by A. Bobylev, Cercignani and G. Toscani in [9], where 
Maxwellian molecules were considered using Fourier transform. A lot of properties of a more 
general model which includes elastic and inelastic collisions for Maxwellian molecules as a 
particular case (but not hard spheres) have also been studied in the chapter of book [5] by A. 
Bobylev, Cercignani and I. Gamba. 

Concerning hard spheres kernel, paper [29] presents a one dimensional model where the 
restitution coefficient depends on the relative velocity of the colliding particles through a 
phenomenological inverse power law. Then, the first important result about tail behaviour 
in any dimension was given by A. Bobylev, I. Gamba, and V. Panferov in [10]. Existence of 
smooth solution to a model with stochastic heating was proven by I. Gamba, V. Panferov, and 
C. Villani in [18J. S. Mischler and C. Mouhot finally proved existence and uniqueness of weak 
solutions to the Cauchy problem for (jl.ip in [26] , The more physical space inhomogeneous 
setting was investigated by D. Benedetto and M. Pulvirenti in [4] for the model introduced 
in [29], but in one dimension of space and velocity. Finally, R. Alonso studied in [1] the Cauchy 
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problem near vacuum for a large class of restitution coefficients, in any dimension of space and 
velocity. 

Let us define more precisely what we call Cauchy problem for (jl.ip . 

Definition 2. Given a nonnegative initial condition fi n G L\c\Q , a nonnegative function / 
denned on [0,T] x W 1 is said to be solution to the Cauchy problem (|l.ip - (|1.2p if 

/€C(0,T;4), 

and if (jl.l|) - (|1.2p holds in the weak sense, namely, 

[ T ( [ f(t,v)dMt,v)dv + (Q e (fJ)(t,-)^))dt=[ f m (v)i>(0,-)dv, 

JO \JR d J JR d 

for any V G C\ ([0,7] x E d ). 

We notice that it is always possible to assume the initial condition in Q, namely of mass 
1 and zero momentum, since we may always reduce to this case by a scaling and translation 
argument, using invariance properties of the collision operator. 

Let us now introduce the cooling time T c of / as the time before blow up, namely 

T c := sup{T > : £ (f)(t) > 0, Vt < T}. 

For a collision kernel given by (jl.4p and f| 1 . 5 [) . Theorem 1.4 of [26] (recalled in the Appendix 
as Theorem IA. lj) states well-posedness in L\ of the Cauchy problem for the granular gases 
equation, together with the existence of the cooling process when t — )■ T c . 

We now state main results of this paper concerning cooling process and blow up solutions to 
the anomalous model. It is known from |26j that if a is greater than 1/2, T c is finite, meaning 
that blow up of / occurs in finite time whereas it takes infinite time for a < 1/2. We give 
the accurate asymptotic behaviour of the energy in the main theorems of this article, for all 
nonnegative values of the parameter a, in agreement with the formal computation (jl.lOp . 

Theorem 1.1. Let the collision kernel be subject to assertions \l-4\ and \ 1.51 If f is solution 
to (jl.ip with an initial condition < fi n G L\ n Q\ H L p for p > 1, there exist some positive 
constants Ci, i G {1,...,6} andT c , depending on e, a and fi n such that 

(i) sub-critical case: if Q < a < 1/2, and then a = l/(2a — 1) < 0, 

(c lt +i)-»» £ £(m 5 ( ci, t +!)-»- Vf =" 0; 

(ii) critical case: if a = 1/2, 

e -C 8 t < £(f)(t) < e~ C4t , > 0; 

(Hi) super-critical case: if a > 1/2, and then a = 1/(2 a — 1) > 0, 

(-C-5 1 + l) 2a < f (/)(*) < (-C 6 1 + l) 2a , Vt < T c . 
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Moreover, for all a > 0, there exist a function V G C 1 (0, T c ) and a nonnegative profile G G 
L\ n Qi such that the distribution 

F(t, v) := V(t) d G(V{t) v), V (t, v) G (0, T c ) x R d , 

is solution to (jl.ip . called self-similar solution. 

We propose in Section [2] a proof of this Theorem based on the analysis of self-similar 
solutions to the inelastic Boltzmann equation with a linear drift term (widely studied in [24, 25J) 
and the asymptotic analysis of the time derivative of the energy. To this end, we introduce a 
new self-similar scaling, nonexplicit and strongly nonlinear (energy-dependant, as in |13j). 



Remark 1 . One can check that the point (i) of this theorem is in good agreement with the 



one proven in |24j with a = 0, and which is classically known as Haff's Law [20]: 



<£(/)(*)< 7—^ 



(l + f)2 -~VA ' ~ (l+t)2- 

We then state a theorem of existence and uniqueness (up to a translation of time) of self- 
similar solutions to equation with a > 0. We also obtain convergence of solutions to (|l.ip 
bounded in L\ toward these particular solutions. 

Theorem 1.2. Let the collision kernel be subjects to assertions (jl.4p and (|1.5p with a > 0. 
There exists a constructive e* G (0, 1) such that for all e G [e*, 1), 

(i) the self-similar profile G = G e G L\ D Q\ is unique and if F e and F e are two self- 
similar solutions to there is a time to < T c such that F e (t,v) = F e (t + to, v) for 
t > max{0, — to}. 

(ii) For p > 1 and any Mo > 0, there exists r\ G (0, 1) such that if 

{ hn g L\ n Q x n lp , 

the solution f to equation (jl.ip satisfies for a nonnegative constant C 

C 



\\f(t,.)-F e (t,-)\\ L i < { 



— if a < 1/2, 
Ce -C if M e t/2 if a = 1/2, 

{C(l-C 6 t) a ^ if a > 1/2, 



where a = 1/(2 a — 1) and /i e = (1 — ry) + 0(1 — e). 

The proof of this Theorem is also based on the analysis of self-similar solutions to the 
inelastic Boltzmann equation with drift term, together with the use of the explicit cooling 
processes of Theorem II. 11 
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1.2 Plan of the Paper 

The article is organized as follows. We prove Theorem 11.11 in Section [2j We give a nonlinear 
rescaling that allows us to treat the three different cooling process at once in Subsection 12. li 
We then use the expression of this rescaling to state a relation between energy in classical 
and self-similar variables. Subsequently, we show in Subsection 12.21 results concerning the 
asymptotic behaviour of solutions to the granular gases equation with drift term, namely that 
their energy and time derivative are uniformly bounded in time. We finally use this result in 
Subsection 12.31 to establish the rigorous cooling process. 

We then apply in Section [3] this precise asymptotic behaviour together with a Theorem of 
uniqueness of self-similar profiles already proven in |25j of S. Mischler and C. Mouhot to show 
the Theorem 11.21 about rate of convergence toward self-similar solutions. 

Finally, we recall in the Appendix some important results we used in the paper. 



This section presents the proof of Theorem 11.11 which studies the cooling process for sub-critical 
(a < 1/2), critical (a = 1/2) and super-critical (a > 1/2) cases. To this end, we will prove 
some new estimates for the granular gases equation, with and without drift term. Let us start 
by introducing some nonclassical self-similar variables. 

2.1 Nonlinear Self-similar Variables 

We shall use a rescaling (seen e.g. in [TB] and |24j ) . in order to prevent the blow up of / by 
"zooming" on the distribution and studying it in self-similar variables. 

We will assume that the collision kernel is given by (|1.4|) and (|1.5|) with a > 0. The granular 
gases equation (jl.ip then reads 



The functions T and V have to be strictly increasing in time for this rescaling to be well 
defined. 

We shall determine T and V (the self-similar variables) and derive the equation followed 
by the distribution g = g(s, w) with s = T(t) and w = V(t) v. The term V(t) d in front of 
g is simply given by mass conservation. Moreover, by making a change of variables in the 
expression of the collision operator and thanks to hard sphere collision kernel, one gets that 



2 Proof of Theorem 11.11 






lim Tit) = lim V(t) = +oo. 

I, t^T c t^T c 



Q e (g(X-),g(X-)) (v) = X'^Q^g^Xv). 
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Especially, the bilinearity of operator Q e together with (|2.2p yields 

(2.3) Q e (f, f)(t, •) = (V(t)) d - l Q e (g, g)(s, ■). 

Remark 2 . Because of (|2.3j) , our approach of the problem cannot be immediately extended 
to the case of variable restitution coefficient. Indeed, if e = e(\u ■ n\) where n is the unit vector 
joining the center of the particles (this case cover the classical visco-elastic hard spheres model, 
see N. Brilliantov and T. Poschel |llj). introducing the self-similar variables f)2.2[) in relation 
(|2~3|) gives 

Qe(f,f)(t,-) = (V(t)) 2 Qe(t,-)(9>9)(s,-), 

that is collision operator becomes time-dependant through a new restitution coefficient e. 
Finally, the new collision equation is not "autonomous" , which prevents us to use some previous 
results on this equation as we will do in the following (e.g. Haff's Law). 

Now, if / follows rescaling ()2.2[) . its time derivative becomes 

(2.4) d t f = (Vityf- 1 (T'(t) V(t) d s g + V'(t) V w ■ (wg)) . 

Thanks to relations (|2.3|) and (|2.4|) . if / is solution to equation (|2.ip . the distribution g is 
solution to 

(2.5) T'{t) V(t) d s g + V'(t)V w ■ (wg) = £(f)(ty a Q e (g,g). 

We shall get rid of the term in (|2.5j) involving a negative power of the kinetic energy to obtain 
a simpler equation, the classical homogeneous Boltzmann equation with an anti-drift term. 

To this end, let us assume the rescaling to be "nonlinear" by asking it to depend on the 
energy of the solution itself: 

V'(t) = r£(f)(t)- a , 
T'(t)V(t) = V\t)/r, 

where r is a nonnegative parameter. This idea is very close to the one used in [13] for nonlinear 
diffusion equations (sometimes referred to as the Toscani Map). It was also extended in [H] 
to show the apparition of chaotic behaviour for a precise (constructive) nonlinearity : the self- 
similar profile of this equation "oscillates" between Gaussian (heat equation) and ZePdovich- 
Kompaneets-Barenblatt (porous medium equation) profiles. 

The functions V and T are not explicit but one can see that they are well defined and 
fulfill the requirement of the scaling (|2,2p . Indeed, the map 1 1— >■ £(f)(t)~ a is strictly increasing 
on [0, T c ) thanks to the cooling process (jA.ip . is 1 when t = given that the energy of the 
initial distribution /j n has been normalized to this value and tends to infinity when t — > T c . 
Moreover, one has T(t) = log(V(i))/r with r > and then T is increasing on [0, T c ) with 
T(0) = and lim t ^ Tc T(t) = +oo. 

With such an expression for V and T plugged in (|2.5p , the distribution g is solution to the 
following equation: 

(2.6) ^+TV w -(wg) = Q e (g,g). 

Actually, this is the granular gases equation for inelastic hard spheres with constant restitution 
coefficient, complemented with an anti- drift term. This last term will act like an input of energy 
and will prevent the blow up of g(s, ■) toward a Dirac mass when s — > oo. This equation has 
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been thoroughly studied in articles \24\ 125] . Especially, we know that the kinetic energy of g 
is bounded by above and below by two positive constants and that there exists steady states 
G for this equation. We summed up this result in the Appendix as Theorem IA.2I 
Summarizing, the distribution / defined by 

' f(t,v) = V(t) d g(T(t),V(t)v), 
< V'(t)=r£(m)- a , 
T{t) = \og{V{t))/r, 



is solution to the granular gases equation (|2.ip for all times as soon as the function g is solution 
to drift-collision equation (|2.6p . Conversely, if / is solution to (|2.ip , one can associate a solution 
g to rescaled equation (|2.6p by setting 

g( S ,w) = e- dTS f{V- 1 (e TS ),e- TS w). 

Note that the inverse V' 1 of V is well defined, by the discussion above. 

Moreover, using the change of variables w = V(t)v, the following relation between energies 
of / and g holds: 

£(g)(T(t)) = (V(t)) 2 £(m). 
Then, the function V can be written for all < t < T c as 

Using this expression, the bounds (|A.2|) of the energy of g and the raw cooling process (|A.ip . 
one has another proof that V(t) —> oo and T(t) = log(V(t)) —> oo when t — > T c . 

Finally, if G is a self-similar profile, that is a stationary solution to (|2.6|) . we may associate 
a self-similar solution F to the original equation (|2.ip by setting 

F{t,v) = V(t) d G{V(t)v), V(t,«) 6 (0,T C ) x R d . 

Such a G exists thanks to Theorem IA.21 which proves the last assertion of Theorem 11.11 

2.2 Preliminary Results 

We will show in this Subsection two results concerning solutions to the granular gases equation 
with and without drift term, that we will need in order to prove rigorously the cooling process. 
Let h = h(t, v) for t > and v G M. d be solution to the inelastic Boltzmann equation 

(2.8) ^ = Qe(h,h), 

that is the collision equation (|2.ip with a = 0. According to [2] or [23], the (sub-critical) Haff's 
Law holds: 

m M 

(2.9) < £{h)(t) < Vt > 0. 

(i + /z *r (i + Mo^r 

An upper control of the third order moment of h can be deduced from this law. Let mi(h) 
denotes the 2Z-th order moment of a distribution h, that is 

mi(h) = / h(v) \v\ 21 dv. 
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Lemma 2.1. Let h be solution to equation ()2.ip for a = ; with an initial condition h{ n £ QC\LP 
for p > 1. There exists a nonnegative constant k such that if m 3 / 2 (hi n ) < k> then for all t > 0, 

(2.10) < m 3/2 {h)(t) < K 

(1 + /J, t) 

Proof. Setting mi = mi(h) and multiplying ()2.ip with a = by \v | 3 , one gets after integration 
in velocity space 

— m 3/2 (t) = / Q e (h,h)(t,v) \v\ 3 dv. 

at J R d 

Some more informations are then needed regarding the third order moment of the collision 
operator. A.V. Bobylev, I.M. Gamba and V.A. Panferov had shown in [10] (and this result 
was then extended to the case of bounded cross section in [2], see Appendix, Theorem IA.4|) 
the following estimate: 

Q e (h, h)(t, v) M 3 dv < -(1 - 7 ) m 2 (t) + 7 S 3/2 (t), 



where < 7 < 1 and 



[|] 

^3/2 = ^2 \ k) ( mfc +V2 m 3/2-fc + m km 3/2 _ k+1/2 ) 



fc=i 

-(m 3/2 my 2 +771?) . 



By definition, we have m\ = £{K). Besides, by convexity, m 2 > Trig* and th\ > m±/ 2 - Thus, 
m 3 / 2 (t) verifies thanks to the sub-critical Haff's Law (|2.9p : 

J t ^s/2(t) < -(1 - 7) m 3/2 (t)^ + ^7 (e (h)(t) 2 + m 3/2 (t) £{h){t) l l 2 

, si/3 3 / M 2 M 1 / 2 



Then, one has 



f rn 3/2 (t)--^)<-(l- 7 )rn 3/2 (t)^ + 3 / '" ' 



*V" 8,4V '' (l + Mf/- V ' ' (1+Mo*) 4 

where k is a nonnegative constant such that 7713/2(0) < k. Then inequality (|2.10p is fulfilled 
for t = and by continuity, the lemma is proved for t < t* with a positive time t*. For t = t*, 
one gets 777.3/2(2*) = «(1 + /Uoi*) -3 , and the last differential inequality reads 



where Ci = 3 (^0 + 7 V M/2J. If the constant k is chosen large enough, the right hand side of 
this inequality is negative, which concludes the proof. □ 
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Thanks to this estimate, we can now show that the time derivative of S(g) is uniformly 
bounded. 

Proposition 2.1. If g is solution to the inelastic Boltzmann equation with drift term (|2.6[) . 
then there exists two constructive constants — oo < Mq < < M\ < +oo such that 

M < ^S(g)(s) < Mt, Vs>0. 
Proof. Given that g is solution to (|2.6p . it is possible to define a new distribution h by setting 

h(t,v) :=Vo(t) d g(T (t),Vo(t)v), 

where we defined as in 



f V (t) = l + /i *, 
\T (t) = log (1 + 

Then h is solution to the inelastic Boltzmann equation (|2.8p . Thanks to the expression of this 
rescaling, one hase 

£(g)(T (t)) = V (t) 2 £(h)(t). 
Differentiating this relation with respect to time gives 

j t £(g) (T (t)) = A(t) + B(t), 

with 

A(t) = 2iM,(l + lMjt)e(h)(t), B(t) = (l + ^t) 2 j t £{h){t). 
Therefore, using the sub-critical Haff 's Law (|2.9p . if t > 0, 

(2.11 — — < i i < — ^ . 

V ; 1 + Mo t - W - 1 + Mo t 

Moreover, by the weak expression (jl.3p of Q(h,h), equation (j2.8|) and ()1 .6j) with a = 0, one 
has on the one hand 

B{t) = (1 + Mo tf / /i) |v| 2 dv 

JR d 

(2.12) = -/3 e (l + Mo0 2 // /i/i*|v-^| 3 dt;(i^ < 0. 
On the other hand, Fubini Theorem and mass conservation yield 

hh* \v — w*| 3 dv dv* < 8m 3 / 2 (h)- 
Then, using inequality ()2.10p of Lemma 12. II in (|2.12p allows to show that B satisfies 

(2.13) _^Pe_ < B , . < 

V ; l + Mo*~ 

Therefore, we can gather inequalities (|2.11l 12.121 12.13P to write 

2.14 VP ° - , - ; < T £ (g) (T t ) < • 

1 + Mo * etc 1 + Mo £ 
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Using the chain rule, we have 

and we can finally write according to (|2.14p . 

8 k B e d _ . , . . 

2m — < —£(g)(s) < 2M. 

fio as 

This concludes the proof. 

□ 

Remark 3 . We notice in the proof that the lower bound Mq is not necessarily a nonnegative 
constant, allowing the energy of g to have some inflection points. 

We are now ready to prove Theorem II .11 
2.3 Anomalous Cooling Process 

We will prove in this Subsection the cooling process of an anomalous gas for a nonnegative 



constant a in the collision kernel (|1.4|) . that is points (i), (ii) and (iii) of Theorem II .11 Both 
upper and lower bounds for the energy of / will be obtained together, thanks to the results of 
Subsection 12.21 

For the sake of simplicity, let us denote for t < T c 

E(t) :=£{f)(t) and E(t) := £(g) (T(t)) , 

and set r = 1. We have already seen in (|2.7p that V(t) = E^) 1 / 2 E(t)~ l l 2 . Then differentiating 
this expression with respect to time yields 



(2.15) 



E{t)~ a = V'(t) 

E'(t) E(t) l / 2 E'{t) 



2{E{t) E(t)) l l 2 2E(t) 3 / 2 ' 
Thanks to the expression of E and a chain rule, one has on the one hand 

E'{t) = ^-(l g{T(t),w)\w\ 2 dw 



dt 

(2.16) =T'(t) X (t), 
where we defined 

X =-£(«(»)) or. 

On the other hand, the time derivative of T is given by 

(2 17) T'(t) - m - 

[IAI) [t) ~ V(t) ~ E{t) l/2 ' 

Finally, gathering relations (|2.15]h (|2.16j) . (|2.17jl and dividing by E(t)~ a , it comes that 

(2.18) - Eit) 1 / 2 E'(t) E{t) a ~^ 2 = 2. 

E(t) 
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Thanks to the result of Proposition 12, 1\ there exists two finite constants Mq < < Mi such 
that Mq < x(t) < M±. Thus, using the uniform bounds (|A.2p of E and the nonpositivity of 
Mq, one has 

Mq < X (t) K 
c ~ E(t) - c ' 

Therefore, gathering equality (|2.18p . bounds (|A.2[> and inequality (jl.lOp imply that 



(2.19) 



where /3 e is given by (|1.7p and £T > because we took Mo < < 2c\. Then, if < a < 1/2, 
the time integration of the two sides of inequality ()2. 19j) between and t and the fact that 
£7(0) = 1 yield 

1 1 

(dt + l)- 2a ~ W - (C7 2 i + 1)- 2 "' 

with a = 1/(2 a - 1) < 0, C\ = -2a K > and C 2 = -2a/3 e > 0. This is the point [(!)] of 
Theorem II .11 Moreover, if a = 1/2, the same argument gives the point (ii) , namely 

exp(-irt) < E(t) < exp(-0 e t). 

Finally, taking a > 1/2 gives the point 



ni 



with C 3 = 2 a K > and C 4 = 2 a P e > 0: 
1-C 3 t) 2a < £(t) < (1 -C 4 t) 2Q - 



Let us now show the result concerning self-similar profiles of equation (jl.ip . 



3 Proof of Theorem 11.2 



For a > 0, we will study in this Section the uniqueness of self-similar profiles of the inelastic 
Boltzmann equation with drift ()2.6|) . The cooling process found in the previous Subsection 
will allow us to state a result concerning convergence of solutions to (jl.ip towards self-similar 
solutions. 

We shall use a Theorem concerning the convergence toward self-similar profiles for small 
inelasticity in the scaled granular gases equation (|2.6p . which has been shown in [25J (we 
recalled it for reader convenience in the Appendix, Theorem IA.3p . For this, let us set 

r = r e := 1 - e, 

in order to balance the dissipation of kinetic energy by the drift. Thanks to this scaling, we 
have uniqueness and attractiveness of the self-similar profiles of equation ()2.6p . 

Let us show thanks to this Theorem the trend to self-similar solution of our problem. We 
have already seen in Subsection 12.11 that if g is solution to (j2.6|) then / is solution to (jl.ip with 

f(t,v) = V(t) d g(T(t),V(t)v), where 



(3.1) 



v(t) _f£( 9 )(T(t))\ 1/2 



T{t) 



log (V(t)) 
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Thus, if G e is the unique self-similar profile of equation (|2.6p . one can find a self-similar solution 
to equation (jl.ip by setting F e (t,v) = V(t) d G e (V(t) v). The uniqueness of this solution up to 
a translation of time can be shown as in [25] to prove the first point of Theorem 11.21 

Moreover, the transformation w —> V(t) v and the rate of convergence toward equilibrium 
(gj]) give 

\\f(t,.)-F e (t,.)\\ Ll =V(t) d [ \g(T(t),V(t)v)-G e (V(t)v)\ dv 

<\\g(T(t),-)-G e \\ L i 
(3.2) < e -(i-^)"«T(*). 

Besides, thanks to the expression (|3.ip of the self-similar variables and the choice r e = 1 — e, 
one has 

v e T{t) = log(F(t)) + 0(l-e). 
The positive lower bound of £ (<7T(t)) i n Theorem lA.2l together with inequality (13 . 2 p and relation 

y ield 

||/(t,-)-F e (t,.)|| L i<C£:(/)(tr/ 2 , 

where yu e = 1 — r] + 0(1 — e) and C is a nonnegative constant. Finally, using the cooling 
process of Theorem 11.11 that we have shown in last Subsection, we can conclude the proof of 
Theorem 11.21 that is the trend to self-similar solution of the solutions to (jl.lj) and the rate 
of convergence depending on the (weak) inelasticity 1 — e and the coefficient a of the energy 
dependent collision kernel (|1.4p . 

Summary and Perspectives 

We have given in this paper the asymptotic behaviour of the space-homogeneous inelastic 
Boltzmann equation for anomalous gases. Depending on a parameter, we can observe in this 
model a blow up in finite time. We quantified the time of blow up and gave the associated 
self-similar profiles. Under a weak inelasticity hypothesis, we also obtained the uniqueness (up 
to a translation of the time) of the self-similar solutions to this equation, and the convergence 
of the classical solutions toward them. 

To prove these theorems, we introduced a new energy-dependant self-similar scaling, which 
leads to the study of the inelastic Boltzmann equation with a linear drift term. We gave some 
results concerning the asymptotic behaviour of the energy of the solutions to this equation, by 
using some well known theorems about this equation. 

Concerning the perspectives of this work, we would like to adapt our results to nonconstant 
restitution coefficient models, such as the viscoelastic one. We also want to prove the rate of 
cooling using a more classical self-similar scaling, generalising the one used in [16] and |24j . and 
moments methods. That would perhaps allow to give up the rather unphysical LP hypothesis 
on the initial condition for the proof of Haff 's Law. 
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A Some Results about the Granular Gases Equation 

In this appendix we will present some important results concerning the inelastic Boltzmann 
equation, with or without drift term, that we used in this paper. 

Let us start by the resolution of the Cauchy problem (jl.ltil.2p for an anomalous gas. 

Theorem A.l ([2B], Theorem 1.4). Let fi n be a nonnegative distribution satisfying 

fin GLlng and f in G BV 4 n l\. 
Then, the following results hold for a cross section given by \1.4\ and \1.5[ 

(i) the cooling time T c is well define and positive; 

(ii) for each T g]0, T c [, there exists a unique solution f G C(0, T;L\) n L°°(0,T;L\) to the 
initial value problem (|l.ip - (|1.2p . Such a solution is nonnegative, mass and momentum 
conservative, and kinetic energy dissipative; 

(Hi) the energy of f is subject to the following asymptotic behaviour 

(A.l) £(/)(*) -> and /(*> ^ $v=o when t -> T c , 

where the convergence of f occurs for weak-* topology o/A4 1 (M rf ). 

We then give a result stating boundedness of solutions to the granular gases equation with 
a drift term, together with the existence of steady-states for this equation. 

Theorem A. 2 ([21], Theorems 1.1 and 1.3). Let gi n G Q n L p for a fixed p > 1 be an initial 
datum for (|2.6|) with constant restitution coefficient e. If g is solution to the associated Cauchy 
problem, then 

(A.2) < c < £{g){s) < ci < oo, Vs > 0. 

Moreover, there exists a self-similar profile < G G L\ n Q : 

TV w -(wG)-Q e (G,G) = 0. 

Then, this result was extended, proving uniqueness and attractiveness of the steady-states, 
under a weak inelasticity assumption. 

Theorem A. 3 (|25j, Theorem 1.1 (i) and (iv)). There exists a constructive e* G (0,1) such 
that for all e G [e*,l), the self-similar profile G e from Theorem \A.S\ is unique and globally 
attractive on bounded subsets of L\: for any M > 0, there exists rj G (0, 1) such that if 

g in G LlnG, \\9in\\L\ < M o, 

then the solution g to equation (|2.6p satisfies 
(A.3) \\g(t,.)-G e \\ Lh <e-^-^\ 
where v e = r e + O (rf ) . 
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The last important result concerns estimates of the moments of the granular collision 
operator Q e . For a nonnegative distribution / and p > 0, we will define the 2p th moment of / 
and Q e {f,f) by 



f(v)\v\*dv, Q p := [ Q e (f,f)(v)\v\^dv. 



This Theorem was first proved for constant cross section and restitution coefficient in [10] and 
then extended to the more general case of Grad's cut-off assumption kernels with variable 
restitution coefficient in [I]. It states that: 

Theorem A. 4 ([10]. Lemma 3 and pQ, Proposition 2.6). For any p > 1, there exists an explicit 
constant 7 = 7 (p, b) 6 (0, 1) such that 

Q P < -(1 - 7)"Vf 1/2 + jS p , 

where S p is given by 



m fp\ 



1/2 m p _ k + m k m p _ k+1/2 ), 



and the binomial coefficients for noninteger p > have been defined by 

p\ _ pip -!)•••(?- k + 1) k>1 /p> _ , 



.kj M ' ~ ' V0 
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